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ABSTRACT 

We present a quantitative spectroscopic study of twenty-seven red supergiants in the Sculptor Galaxy 
NGC 300. J-band spectra were obtained using KMOS on the VLT and studied with state of the art 
synthetic spectra including NLTE corrections for the strongest diagnostic lines. We report a central 
metallicity of [Z]= —0.03 ± 0.05 with a gradient of —0.083 ± 0.014 [dex/kpc“*], in agreement with 
previous studies of blue supergiants and H ii-region auroral line measurements. This result marks the 
first application of the J-band spectroscopic method to a population of individual red supergiant stars 
beyond the Local Group of galaxies and reveals the great potential of this technique. 


1. INTRODUCTION 

The chemical enrichment of galaxies with heavy ele¬ 
ments lays a complex roadway across the age of the star¬ 
forming universe, and techniques to measure metallicity 
are the astronomer’s compass and sextant as the evo¬ 
lution of galaxies is mapped. The primary tools of the 
trade so far have been the studies of the strongest emis¬ 
sion lines of H ii regions. They have provided impor¬ 
tant first hints at the enrichment history of the observ¬ 
able universe which include gradients across the disks of 
galaxies, a relationship between galactic mass and cen¬ 
tral metallicity, and metallicity evolution as a function of 
redshift. Still, these strong-line H ii methods are plagued 
with systematics arising from their empirical calibrations 
and the complexity of the systems p roducing the effi¬ 
ciently observable emission lines used dKewley & : Ellison 
2008 Kudritzki et al.||2008 Bresolin et al. 12009). At the 
same time, it is a basic tact that the stellar population 
of galaxies provides the drive for chemical enrichment 
as well as the radiated light used to measure that pro¬ 
cess. Thus, the natural markers of metallicity are then 
the stars themselves. 

Only recently have developments in the modeling of 
stellar atmospheres, observational techniques, and statis¬ 
tical methods added the quantitative spectroscopy of su¬ 
pergiant stars to the extragalactic astronomer’s toolset. 
This pioneering work provided a measurement of the 
radial abundance gradient of NGC 300 using optical 
spectroscopy of blue supergiant stars and grids of syn¬ 
thetic spectra calculated without assuming local ther- 
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modynamic equilibrium (LTE). This non-LTE (NLTE) 
quantitative technique revealed a central metallicity of 
slightly below solar and a c lear gradient across th e star¬ 


forming disk of the galaxy (Kudritzki et al. 


2008). 


become 


Since that early work, blue supergiants have 
a powerful tool for measuring metallicities, gradients. 
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chemical abundances obtained from the nebular emission 
line method based on the determin ation of the gas tem - 
perature via the weak auroral lines dBresolin et al.|2009 ) 


albeit most recent work (Kudritzki et al.||2013 ' 2014) in¬ 


dicates a small systematic difference of ~0.1 dex. By 
comparison, metallicities derived from any of the host of 
strong line H ii methods return, at best, consistent trends 
with huge offsets in overall enrichment level. 

A growing body of observational work demonstrates 
that the quantitative spectroscopy of red supergiant 
(RSG) stars rivals the precision of metallicity measure¬ 
ments using blue supergiants and is applicable over sim¬ 
ilar distance scales using existing telescopes and in- 
struments. The pi lot near-infrared J-band study of 
Davies et al. ([201 ^ and hig h spectral resolution fol- 


lowup by Cazak et al. (2014) —both targeting galactic 
RSGs—demonstrate the applicability of the techn ique to 
stars with ro ughly solar chemical enrich ment. In |Davi^ 


et al. (2014) and |Patrick et al. (2014), the technique 
IS successfully applied to RSGs signihcantly below solar 
chemical abundance in the Small and Large Magellanic 
Clouds and the metal poor Local Group dwarf galaxy 
NGC 6822. 

These studies are aided by continued development of 
strongly improved quantitative synthetic spectra used to 
extract stellar parameters. While the underlyi ng grid of 


MARCS stellar atmospheres is calculated in LTE (Gustafs- 


son et al. 2008), corrections have been developed tor 


INL'l'E line formatio n of the strongest diagn ostic lines 
(iron and titanium — ^ergemann et al. (2012 ), silicon — 


Bergemann et al. (20 
et al.M2014| ). -^ 


^ ), and magnesium — Bergemann 
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The RSG technique will become more important with 
the next generation telescopes such as the TMT and 
E—ELT. These telescopes and their instruments will be 
optimized for observations at infrared wavelengths, using 
adaptive optics supported multi object spectrographs. In 
consequence, investigating stars radiating strongly in the 
near-IR—including red giant, asymptotic giant branch, 
and red supergiant stars —will have a clear advantage in 
the future. 

In this paper we report another significant step. We 
have used the multi-IFU NIR spectrograph KMOS at the 
ESO VLT to test the J-band RSG technique over a wide 
range of chemical enrichment across the sta r forming disk 
of NGC 300 at a distance of 1.88 Mpc (Gieren et al. 


20051. 

This first application of the technique beyond the Lo¬ 
cal Group of galaxies is yet another milestone towards a 
large scale application of the RSG technique to the un¬ 
derstanding of the present day chemical enrichment of 
galaxies undergoing star formation. 

After discussing the KMOS observations in we out¬ 
line our analysis procedure in ^ Eor a ful l description 
of the J-band technique we direct readers to |Gazak et al.| 
(2014). Results and discussion are presented in ^ 


2. OBSERVATIONS 
2.1. Target Selection 

We prepared an initial database of RSG candidates us- 
ing the AGS Nearby G alaxy Survey Treasury (ANGST: 
Dalcanton et al.||2009), a public database of stellar pho- 
tometry obtained with Hubble Space Telescope obser¬ 
vations. The ANGST catalog contains 6 fields in NGC 
300 from the Advanced Camera for Surveys (ACS) and 3 
fields observed with the Wide Field and Planetary Cam¬ 
era 2 (WFPC2). The layout of these fields can be seen 
in Figure We select the red F814W (“I”) filter for 
magnitude cuts. The color cut for both instruments 
was “V-I”, using F555W and F814W filters for ACS, 
F606W and F814W filters for WFPC2. The cool temper¬ 
atures and extreme luminosities of RSGs allow for high 
fidelity selection using such color-magnitude cuts. This 
is demonstrated in Figure the “RSG plume” (boxed) 
separates from the population of fainter red objects in m/ 
magnitude dimension and is distinct from other bright, 
hotter objects with a red V—I color. Using the ANGST 
photometry we define RSG candidates as those objects 
with colors bounded in V—I color between 1.2 and 3.5. 
and having m/ brighter than 19.5, a value selected to 
achieve our tyget signal to noise ratio (SNR) of ~100 
(Gazak et al. 2014). Seeking the best spectra possible 
we divided these targets into two lists separated at m/ 
> 19.0 and gave higher priority to the brighter group. 

We note that candidate objects within the fainter 
group could be, in principle. Asymptotic Giant Branch 
stars (AGB). Stellar evolution models predict some over¬ 
lap between RGB and so-called super-AGB at log L/Lq 
^ 4.5. However, these objects would also be young with 
ages smaller than 100 Myr with abundances of Fe, Si, 
Ti unaltered from their initial composition so that they 
would also be good tracers of the host galaxys present 
day abundances. 

Because the HST coverage of NGC 300 is incomplete 
we used overlaps in our database with a more com¬ 


plete but shallower cat alog of B and V photometry 
(Pietrzyhski et al. 2001) to train a B—V vs my color- 
magnitude cut. Candidates selected by this method were 
placed in a third ranked priority group such that HST 
selected targets would be observed first but that no IFUs 
would go unused. 

As a final adjustment to our RSG candidate list we 
search for overlaps between our exist ing candida t es an d 
Spitzer photometry of NGC 300 from Khan et al. (|2010 ). 
The important Spitzer diagnostic for RSG candidacy .s 
the color-magnitude plane defined using IRAC Band 1 
(m 3 6) and the color Band 1 — Band 2 ([3.6] —[4.5]). 
RSGs separate as the brighte st m.^.^ targets with color 
blue ward of 0.0 i n [3.6] —[4.5] (Verhoelst et al.|2009| Bo-|| 
nanos et liripOOO ). We upgrade ail objects in our candi¬ 
date list with overlap as Spitzer RSG candidates to the 
top priority group. Coordinates and photometric data of 
our targets finally used for the analysis are provided in 
Table [5j 

2.2. Observations 

The resulting candidate catalog was used as input to 
the KMOS ARM Allocator (KARMA) software for plan¬ 
ning th e setup of the 24 KM OS Integral Field Units 
(IFUs) (jWegner & Muschielok|]2008j) . These 2".8 square 
IFUs can be placed over the 7'.2 KMOS field of view 
(shown in Figure [^. We observed an inner and outer 
field to cover a significant range of the radial extent of 
NGC300. These fields were observed in 600s integrations 
over the nights of 2013 October 14, 15, and 16. The 14th 
and 15th nights were clear and stable at Paranal with 
median seeing of «l".l. Weather degraded slightly on 
the night of the 16th with intermittent cloud cover and 
variable seeing with the median value around 1".3. 

KMOS was operated in nod to sky mode with science 
integrations of 600 seconds. Telluric standards were ob¬ 
served down each IFU at a frequency of once per every 60 
to 90 minutes. KMOS IFU data cubes were flat fielded, 
wavelength calibrated, and telluric corrected using the 
KMOS pipeline (vl.3.2) provided publicly as a specific 
instance of the European Southern Observat ory (ESO) 


Reflex auto mated data reduction environment (Freudling 

eraI][M^ . 


We find inhomogeneities in spectral resolution across 
the field of view of each KMOS IFU as shown in Fig¬ 
ure [Sj To correct for this we first measure the effective 
spectral resolution of night sky lines in each spatial pixel 
of the inner 12 x 12 IFU field. Each sky line is well ap¬ 
proximated by a gaussian and for each spatial pixel we 
calculate the weighted mean of the full width half max¬ 
ima of those gaussian models. The top panel of Figure]^ 
shows one such example. The effective resolutions, in 
this case, vary by up to AR~300. We correct this by 
smoothing the sky and science spectra for each IFU to 
bring the field to a more constant resolution of AR~100 
lower than the minimum measured effective resolution. 
The results of re-measuring effective resolution are dis¬ 
played in the bottom panel of Figurej^ which has a global 
offset in measured resolution but is jotted on the same 
dynamic scale. The effective resolution after these cor¬ 
rections ranges from the expected value after smoothing 
to values smaller by AR^IOO. 

Once the data cubes are homogenized we extract a sci¬ 
ence spectrum. We do this by summing over each wave- 
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Fig. 1. — The spiral galaxy NGC300. The Hubble ACS fields are overplotted in black, with WFPC2 frames in black dashed with gray. 
Observed RSGs are circled in red. In the upper left a gray circle denotes the size of the KMOS field of view, within which 24 IFUs can be 
placed per pointing. 
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Fig. 2.— Color Magnitude selection for extragalactic RSGs. The 
left panel shows color vs magnitude for the WFPC observations of 
NGC 300 while the right panel is the same but for ACS fields. The 
boxed regions show the brighter “priority one” RSG candidates in 
this color magnitude space and the fainter “priority two” region 
below. 

length slice in the IFU data cube. In practice, we find 
that the highest signal to noise is recovered by extract¬ 
ing in a two pixel radius around the peak target flux in 
the science frame. We suspect that inhomogeneities in 
wavelength across the held of view play some role in this. 



After constructing ID spectra of object and sky frames, 
we utilize the ESO tool Skycorr which takes as input 
those two spectra and scales airglow lines using theoret¬ 
ical knowledge of the night sky lines. The resulting sky 
spectrum provides a removal of variability in sky emis¬ 
sion based on diff erence in time a nd sky position of the 
two observations (Noll et al.||2014 ). 

As each held was observed multiple times, the hnal 
science spectra were produced by extracting the median 
of all individual sky-subtracted spectra. This spectral 
atlas is plotted in Figure]^ 


2.3. Target validation 

In order to verify that our selected targets were in¬ 
deed members of NGC 300 and not foreground objects, 
we studied the radial velocities Urad of each of the stars. 
We measured Urad by cross-correlating the reduced spec¬ 
tra with a model spectrum from our grid. The choice of 
model spectrum was found to be unimportant, produc¬ 
ing variations in the measured Uj-ad of less than 1km s“^, 
insignihcant compared to the other sources of error (see 
below). The Urad were then corrected for barycentric mo¬ 
tion and converted to the Local Standard of Rest (LSR) 
using the Starlink package RV. 

For the errors on Vmd we hrst specihed a lower limit to 
the uncertainty of 1/lOth of a pixel, which at the 2-pixel 
sampling of KMOS corresponds to ±5kms“^. To this 
error we added in quadrature the observed spatial vari¬ 
ations in wavelength calibration across each IFU, which 
we found to be ±6kms“^. This gives us an error on each 

Urad of ±8kmS“^. 

To determine what the expected Urad would be 
position of each target, we used the Hi map of 
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Fig. 3. — Variation of spectral resolution of the inner 12x12 spa¬ 
tial pixels of a typical KMOS IFU measured from night sky lines. 
The top panel shows a typical measurement after a 600 sec expo¬ 
sure. The bottom panel presents the resolution after the homoge¬ 
nizing procedure described in the text is applied. 


meier et al. (|2011), and their model of NGC 300’s rota¬ 
tion curve, 'i'hese authors fit the observed velocity map 
with a model which allowed the rotational velocity Urot, 
inclination i and position angle PA to vary with galac- 
tocentric distance. Since the variation in i was found to 
be small (45±5°), we kept this parameter fixed at 45°. 
For Urot (=40-100kms“^) and P.A(=290-330°), we used 
fits to their observed trends with galactocentric distance. 
The predicted radial velocity Upred as a function of posi¬ 
tion on the plane of the sky relative to the centre of the 
galaxy (x, y) was then determined from, 


'^pred(4^5 2/) — 


'sys 


- Urot COS 0 sin i 


( 1 ) 


where Vgys = 136km s 
cos 9 is given by, 


-1 


(Westmeier et al. 20111, and 


and. 


( —) 
Vcosz/ 


( 2 ) 


(3) 


In Fig. we plot the observed versus predicted ra¬ 
dial velocities for each star in our sample. The Hrst 
thing we can say is that there are no obvious candidates 
for foreground stars in our sample, since such objects 
would be expected to have low radial velocities of a few 
xl0kms“^. For the majority of our sample, there is ex¬ 
cellent agreement between their measured radial veloci¬ 
ties and the model rotation curve to within la. The only 
deviations greater than 2tT are seen at high Urad, where 
four datapoints lie ~20kms“^ above the predicted trend. 
If NGC 300 is similar to the Milky Way in that devia¬ 
tions from the model rotation curve of 20km s“^ are not 


uncommon (|Russeil [2003|) , then we can say all stars in 
raOL 


our sample 


NGC 300 members. 


lave radial velocities consistent with being 


3. ANALYSIS 


Observed spectra are analyzed through a comparison 
with synthetic spectr a based on a grid of ID L TE MARCS 
model atmospheres (Gustafsson et al. 20081. The line 
formation calculations of the synthetic spectra are car¬ 
ried out in NLTE for the elements iron, titanium, silicon 
and magnesium, which produce the most critical diagnos¬ 
tic lines prominent in the J-band spectra. The details of 
the NLTE radiative transfer method and atomic models 


are described in Bergemann et al. (2012,2013 2014). All 
other lines including the (weak) molecular contributions 
are calculated in LTE. The free parameters of this grid 
include effective temperature, log gravity, metallicity [Z] 
(normalized to Solar values, [Z] = log Z/^), and micro¬ 
turbulence (Teff, log 5 , [Z], ^ — see Table [Tfor parameter 
coverage). For the analy sis we apply the sam e technique 
as described in detail in |Gazak et al.| ( [2014[ ). We note, 
however, that we use an updated set ot synthetic spec¬ 
tra calculated with i mproved oscillator strength for the 
J-band Til lines (see Davies et al. 2014). 

For a discussion of the limitations and possible im¬ 
provements of the model atmospheres and the radiative 
transfer methods used we refer to the publications cited 
above and the references therein. On the other hand, 
after the successful tests of the method at Milky Way 
( GnzaketalJ|201d| and significantly l ower metallicities 
(Davies et al.||2014| |Patrick et al.|2014) we are confident 
that the method drivers reliable metallicities. 

In addition to model parameters we fit for “effective 
resolution”, Rqs, an unresolvable combination of the ef¬ 
fects of instrument spectral resolution R (JA/A) of each 
IFU and a very small contribution of macro turbulence 
in the atmospheres of RSGs. This is done by degrading 
each model to a set of resolutions ranging from twice the 
expected spectral resolution downwards until a clear 
minimum has been defined. The minimum of a parabolic 
fit through the vs. R is adopted as the best ht R^s 
for that combination of model and data. After an ini¬ 
tial construction of the x^ grid allowing each model to 
be measured with its own ideal i?off we lock the i?off to 
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Fig. 4.— NGC 300 RSG spectra plotted in black with corresponding best model fits in red. Diagnostic features are marked. Each object 
; labeled to the right of the plotted axis and information corresponding to the fit is tabulated by this name in Table 
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Fig. 5. — The observed radial velocities of our sample stars versus 
that predicted u s ing th e model rotation curve of NGC 300 from 
|Westmeier et al!] ||2011[l. 


TABLE 1 

MARCS Model Grid 


Parameter 

Notation 

Min 

Max 

Spacing 

Eff. Temperature [K] 

Teff 

3400 

4000 

100 



4000 

4400 

200 

Log gravity 

logs 

-1.0 

-tl.O 

0.5 

Metallicity [dex] 

[Z] 

-1.00 

-1-1.00 

0.25 

Microturbulence [km/s] 


1.0 

6.0 

1.0 


Note. — Parameter grid for marcs atmospheres (and 
synthetic spectra) utilized in this work. 


We measure best fit parameters by isolating the six 2D 
planes in each parameter pair combination which contain 
the minimum As a result, two parameters are locked 
to the best fit values in each of the six slices. We inter¬ 
polate the grid of each plane to a parameter grid with 
four times the density and take the minimum of that 
dense grid as the best fit values for the two parameters 
defining the plane. After completing the procedure for 
all six planes, the three measurements of each parameter 
are ave raged into a set of best fit parameters (again we 
refer to|Gazak et al.||2014|for all details). 

Pararneter uncertainties are assessed with a Monte 
Carlo simulation. This is accomplished by adding 1000 
instances of random gaussian noise to a linear interpo¬ 
lation from the model grid to the best fit parameters. 
The gaussian noise is scaled such that the SNR of this 
noisy model matches the SNR of the observed spectrum. 
For each instance we determine fit parameters and for 
each parameter we define a Icr uncertainty such that the 
central 68% of the 1000 measured parameters lies within 
±1(7. The best fit parameters and their derived uncer¬ 
tainties are presented in Table 


4. RESULTS AND DISCUSSION 
4.1. Results 

The parameter fits are in general quite precise with 
individual metallicities as accurate as 0.10 to 0.15 dex. 
The exception is object 134 with a high effective temper¬ 
ature just at the edge of our grid of synthetic spectra. 
It may well be that this object is slightly hotter which 
would reduce its metallicity and improve the fit of the 
silicon lines. 

In the following, we discuss metallicity and metallicity 
gradient and the evolutionary status of our objects. The 
main results are summarized in Figures 


the value which matches the minimum of that 
re-calculate. This step is completed for accurate mea¬ 
surement of parameter uncertainties. 

Correcting the continuum level of our models to that 
of the observed spectra starts with selecting the high¬ 
est model flux point in wavelength bins of characteristic 
“continuum width” CW = X/2Reff such that the fit 
samples the full spectrum. We note that at the low reso¬ 
lution of spectra in this work, the true continuum is lost 
to a blended forest of molecular lines such that we are 
actually fitting a “pseudocontinuum”. While the level 
of the pseudocontinuum is certainly a function of the 
physical stellar parameters, especially [Z], we find that 
no detrimental effects due to this phenomenon appear to 
resoluti ons significantly bel ow that of our observed spec¬ 
tra (see Cazak et al. 20141. After removing continuum 
outliers we fit a smooth third order polynomial through 
the ratio of observed flux to model flux at the continuum 
points as a function of wavelength. This function is then 
applied to scale model to data. 

With the continuum normalization procedure de¬ 
scribed in the preceding paragraph it is important to 
realize that our spectroscopic J-band method is not af¬ 
fected by the potential presence of circumstellar dust. 
Though such dust may add a small amount of extinction, 
it will affect the pseudo continuum and the lines equally 
and, thus, leave the abundance measurements unaltered. 


4.2. Metallicity Gradient 

This first measurement of metallicity of the young stel¬ 
lar population of stars across the disk of a spiral galaxy 
beyond the Local Croup using the RSG J-band method 
requires a careful comparison with previous work. 

The most direct comparison of our met allicity results 
is to the work by Kudritzki et al. (2008), who studied 
24 A- and B-type supergiants across the radial extent 
of this galaxy. They used synthetic NLTE spectra to 
derive metallicities from an ensemble of heavy atomic 
species (Mg, Si, S, Ti, Cr, and Fe) and derived logarith¬ 
mic metallicities relative to the Sun. We display their 
results in Figure ]C With a linear regression we obtain 
a central metallicity of —0.07 ± 0.09 dex and a gradient 
of —0.081 ± 0.011 dex kpc ~^ from this samp le of BSGs. 
We note that the original Kudritzki et al. (2008) val¬ 
ues for central metallicity and gradient were —0.06 dex 
and —0.083 dex kpc“^, respect ively. We applied the n ew 
galactic orientation model of Bresolin et al. (2009|) to 
calculate slightly different galactocentric distances. This 
resulted in a small change in the central metallicity and 
metallicity gradient. _ 

Our J-band met hod is a natural extension of the IKu-l 
dritzki et al. (2008) technique and the results for our ob- 
served RSC sample are also shown in Figure!^ We stress 
that we use an independent but evolutionarily connected 
population of stars, as blue supergiants below masses of 
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Fig. 6. — The stellar metallicities and metallicity gradients in NGC 300 at 1.88 Mpc. Blue squares are BSGs from |Kudritzki et al. ]( |2008| l 
and a linear regression through those points is plotted as a blue dashed line. Black circles and the corresponding solid black linear regression 
line are RSGs from this work. 


TABLE 2 

NLTE Stellar Parameters eor NGC 300 Red Supergiants 


Target 

R/R 25 

Teff [K] 

logs 


Z [dex] 

? [km 

-/s] 


010 

0.069 

4096 

± 

40 

-0.88 

± 

0.2 

- 1 - 0.01 

± 

0.09 

4.4 

± 

0.4 

3500 

022 

0.120 

3535 

± 

70 

-0.73 

± 

0.3 

-0.14 

± 

0.11 

3.1 

± 

0.2 

2800 

024 

0.204 

4153 

± 

40 

-0.16 

± 

0.2 

-0.16 

± 

0.12 

4.0 

± 

0.2 

3000 

007 

0.250 

3912 

± 

70 

-0.52 

± 

0.2 

-0.12 

± 

0.09 

4.0 

± 

0.2 

2700 

013 

0.282 

3875 

± 

50 

-0.47 

± 

0.5 

-0.22 

± 

0.15 

4.3 

± 

0.2 

3500 

127 

0.293 

4220 

± 

30 

-0.80 

± 

0.3 

-0.10 

± 

0.10 

4.3 

± 

0.3 

2800 

Oil 

0.330 

4063 

± 

70 

-1-0.47 

± 

0.2 

-0.09 

± 

0.15 

4.1 

± 

0.2 

2800 

009 

0.332 

4075 

± 

40 

-0.12 

± 

0.2 

-0.18 

± 

0.12 

4.3 

± 

0.2 

3100 

014 

0.384 

4260 

± 

30 

-0.38 

± 

0.2 

-0.25 

± 

0.08 

3.6 

± 

0.2 

3400 

006 

0.432 

4205 

± 

90 

-0.42 

± 

0.3 

-0.31 

± 

0.10 

4.1 

± 

0.2 

3200 

101 

0.487 

4228 

± 

42 

-0.52 

± 

0.1 

-0.22 

± 

0.09 

4.2 

± 

0.2 

2800 

102 

0.551 

4010 

± 

147 

-1-0.07 

± 

0.2 

-0.31 

± 

0.10 

3.9 

± 

0.2 

2800 

133 

0.564 

3953 

± 

52 

-0.40 

± 

0.1 

-0.28 

± 

0.12 

3.6 

± 

0.2 

3500 

130 

0.569 

4123 

± 

30 

-0.44 

± 

0.2 

-0.24 

± 

0.11 

4.1 

± 

0.2 

3000 

126 

0.589 

4251 

± 

60 

-0.41 

± 

0.2 

-0.21 

± 

0.12 

4.1 

± 

0.2 

2900 

135 

0.590 

3953 

± 

53 

-0.32 

± 

0.2 

-0.30 

± 

0.11 

3.8 

± 

0.3 

3200 

134 

0.610 

4400 

± 

20 

-0.17 

± 

0.7 

-0.15 

± 

0.29 

5.0 

± 

0.2 

3400 

039 

0.647 

4240 

± 

50 

-1-0.56 

± 

0.3 

-0.42 

± 

0.11 

3.0 

± 

0.2 

3000 

128 

0.672 

3776 

± 

63 

-0.53 

± 

0.4 

-0.17 

± 

0.12 

5.0 

± 

0.2 

3000 

129 

0.776 

4113 

± 

50 

-0.28 

± 

0.2 

-0.45 

± 

0.09 

4.1 

± 

0.3 

2900 

139 

0.866 

4280 

± 

30 

-1-0.56 

± 

0.4 

-0.42 

± 

0.10 

4.1 

± 

0.2 

3100 

034 

0.918 

3833 

± 

44 

-0.01 

± 

0.2 

-0.48 

± 

0.09 

4.1 

± 

0.2 

2500 

055 

0.949 

3981 

± 

119 

-0.40 

± 

0.1 

-0.45 

± 

0.17 

4.1 

± 

0.2 

3100 

132 

0.962 

4160 

± 

43 

-0.26 

± 

0.2 

-0.34 

± 

0.12 

4.5 

± 

0.2 

2800 

140 

0.970 

4141 

± 

30 

-1-0.04 

± 

0.3 

-0.42 

± 

0.11 

4.1 

± 

0.2 

3200 

031 

0.975 

4270 

± 

40 

-1-0.20 

± 

0.2 

-0.41 

± 

0.12 

4.0 

± 

0.2 

3400 

037 

1.014 

3912 

± 

50 

-1-0.40 

± 

0.2 

-0.53 

± 

0.14 

3.1 

± 

0.3 

2900 


Note. — R 25 — 5.33 kpc (Bresolin et al., 2009) based on a distance of l.SS Mpc JGieren| 
|et al |2005[l. 


CO 3 OM 0 evolve into the RSGs we observe as they exhaust 
their core hydrogen. By all means, then, the metallici¬ 
ties obtained with the two techniques should agree. Still, 
the stellar models and synthetic spectral calculations are 
independent for our two techniques as the physical con¬ 
ditions of these two populations are distinct. The cool, 
inflated atmospheres of RSGs contain a host of singly 
ionized metals and a forest of molecular features which 
are absent in the hot spectra of BSGs. We also observe 


our stars in a different wavelength regime subject to sep¬ 
arate observational difficulties. Finally, due to the added 
complexity of RSG atmospheres it is not yet possible to 
perform calculations in pure NLTE, even though our grid 
of synthetic spectra do have NLTE calculations for the 
most critical spectral features. Despite these differences, 
the J-band technique derives metallicity as a ratio of the 
solar abundance pattern just as the BSG method does, 
using strong atomic lines of Fe, Ti, Si, and Mg. Our 
























































observations of NGC300 RSGs then become an excellent 
test of both methods and of the theoretical calculations 
from which they draw. 

Using our observations and the RSG J-band technique, 
we measure a central metallicity of —0.03±0.05 dex and a 
gradient of —0.083±0.014 dex kpc“^. This is a stunning 
agreement between BSGs and the RSG J-band method 
a nd certainly repr e sents a breakthrough of the technique. 


Bresolin et al. (20091 provide an important second 
measurement of Ihe central metallicity and gradient of 
NGG 300 using measurements of Hii region auroral lines. 
Studies using this “direct method” of measurement agree 
with abundances derived from populations of stars in en¬ 
vironments wh ere metallicity is below solar, for example, 
dwar f galax ies ( Stasihsl^|2005 Bresolin et al.|[2006 Lee 


et al. 


2006), and for some regions in the Milky Way and 


M33 (MW: Rolleston et al. 2000 Deharveng et al. 2000 

Daflon & Gunhal 

2U041 M33: IV Ichez et al.| 

19881 |Ur-| 

baneja et al.||2005 

Bresolin||2011 ). In this tec. 

inique the 
;al param- 
isentangle 
ibundance 
a central 
3.02 and a 

[UiiiJ aurora 
eter, the elec 
the effect of 
and tempera 
oxygen abun 

1 lines give access to a key physic 
tron temperature T^,, which can c 
line strengths based on oxygen 

ture. 

Bresolin et al. (2009) fine 

dance ot 12-kiog(0/H) ^^7 ± 


gradient of —0.077 ± 0.006 dex kpc“^. The gradient is 
slightly shallower than for the RSGs and BSGs, but the 
difference is small. The comparison of central metallic¬ 
ity depends on the assumed value for the sola r oxygen 
abundance . Choosing 12-|-log(O/H)0 = 8.69 (Asplund 
et al.|2009 ) returns a H ii region central metallicity value 
ot Z = —0.12 ± 0.02. While this value agrees with the 
BSG result it indicates a small offset (~0.09 dex) relative 
to the RSGs. We note, however, that such an offset is 


also f ound in the recent work by Kudritzki et al)] (|2013 
20141 in the spiral galaxies NGG 4258 and 3621 wheii 


BSG and Hii region metallicities are compared. A 
possible reason for this discrepancy is the depletion of 
ne bular oxygen on dust gra ins as, for instance, suggested 
by Zurita & Bresolin (20121. On the other hand, we note 
that comparing the metallicities obtained with the three 
methods at R/R 25 = 1 the situation looks much better. 
We find [Z] = -0.50, -0.47, -0.53 for the BSG, RSG, 
and H ii auroral line methods, respectively. Given that 
the three methods are entirely different the agreement is 
striking. The measurements from this work are tabulated 
with the discussed work from the literature in Table [3| 

4.3. Stellar Evolution 

Stellar evolution theory predicts a shift of the Hayashi 
limit (the minimum temperature of an RSG evolution¬ 
ary track, when stars are nearly fully convective) towards 
higher effective temperature with decreasing metallicity. 
For a stellar evolutionary track with solar metallicity 
and a mass of M = 1 5 MfT> a minimum temperature 


of ~3600K is obtained (Ekstrom et al. 20121 Meynet & 
Maeder 2002[), whereas models with SiVlG like metallic¬ 


ity (jGeorgy et al.|2013 ) have a minimum temperature of 


41UUK. In contrast to this prediction, the RSG J-band 


6822 

Gazaic et al.|2014 Davies et al.|2014 

Patrick et al. 

2014) 

does not show such a trend and leac 

s to the con- 


elusion that the temperatures of RSGs are independent 
of metallicity. We can use our observations in NGG 300 
as an additional test in this regard since we encounter 


a range of metallicities between solar in the center of 
the NGG 300 disk down to 173 solar at the isophotal 
radius. Accordingly, Figure 0 shows a plot of effective 
temperature against metallicity. It confirms the absence 
of a trend at least in the metallicity range covered. In 1- 
D evolutionary models, the temperature of the Hayashi 
limit (i.e. the maximum size of a star) is governed by 
the treatment of convection, in particular by the mixing- 
length parameter a. The value of a is typically tuned in 
order to reproduce the properties of the Sun, but then is 
fixed at this value for stars of all masses, ages and metal¬ 
licities. This fixed value of a typically results in evolu¬ 
tionary paths for RSGs which are ‘inclined’, decreasing 
in Teff as L increases, as well as producing trends of lower 
Teff for RSGs with higher initia l masses, and higher av- 
erage T^fp at lower metallicities (iMeynet & Maeder|2000 

-I—1 • 1 . 1 II /~V I I -r-> , , . 1 *1IT ^ J 7^ ^ H • 


Eldridge et al.|[2008 Brott et arj|2Ull ). 


Recently, studies using 3-U hydrodynamical simula¬ 
tions of convection have attempted to empirically model 
the evolution of a for stars across a. range of TpfF, log q 
and m etallicity (Trampedach et al. 2014 Magic et al. 


20151. The nature of the simulations used in these stud- 


les meant that it was not possible to cover the parameter 
space occupied by RSGs, and instead focussed on stars of 
higher gravities and temperatur es. However, extrap olat- 
ing from the observed trends in Magic et al. (2015), the 
results seem to suggest that a should increase at lower 
metallicity for a fixed Tea and log g. If correct, this would 
cause RSGs at a fixed initial mass and evolutionary stage 
to become hotter at lower metallicity, exacerbating the 
trend seen in evolutionary models rather than counter¬ 
acting it. In effect this is the opposite to what we find 
in our results, that the average effective temperature of 
RSGs seems to be independent of metallicity. 

The discussion of Figure is, of course, based on the 
assumption that our determination of effective tempera¬ 
tures using the J-band technique is not affected by spu¬ 
rious systematic effects which lead to erroneous effective 
temperatures which are always very similar even for sig¬ 
nificantly different metallicities. A first clear in dication 
that this is not the case comes from the work by |Davi^ 


et al. (2013). This comprehensive study of RSGs in the 
LMG and SMC determined effective temperatures us¬ 
ing ESO VLT XShooter observations of SEDs covering 
the full spectral range from the near-UV to the near- 
IR. It was found that the RSG effective temperatures in 
both galaxies are very similar and do not seem t o be af¬ 
fected by th e difference in metallicity. Moreover, |Davies| 


et al. (2014) recently re-investigated the same sample of 
RSGs applying the spectroscopic J-band technique and 
found effective temperatures consistent with the fits of 
full SEDs indicating that the J-band method provides 
robust estimates of effective temperatures. For a more 
detailed di scussion of RSG tem per ature determination s 
we refer to Davies et al. (2013) and Gazak et al. (2014). 

We construct a Hertzsprung-Russell diagram (HRD) 
to compare our results with stellar evolution theory. We 
calculate bolometric luminosities for program stars us¬ 
ing F814W (I) band photomet ry from Hubble an d the 
bolometric correction recipes of Davies et al. (2013). We 
adopt t he distance modu l us of m = 26.37 ±1105 trom th e 
work by Kudritzki et al. (2008) and Gieren et al.j (|2005|). 
These luminosities are plotted against the ettective tem¬ 
peratures from our spectral fit in the HRD of Figure 
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Fig. 7.— Effective temperature of our RSG targets as a function 
of metallicity. 



Ta,[K] 


Fig. 8 .— Hertzsprung-Russell diagram of program stars for which 
we have F814W (I) band photometry from Hubble. Temperatures 
are the parameter fits from our analysis procedure and we calculate 
luminosities using the bolometric correction recipes for RSGs of 
[Davies et al.| |[2013|l. 


We then overplot evolutionary tracks with solar metal¬ 
licity adopting the Geneva database of stellar evolution¬ 


ary models in cluding the effects of rotation (Meynet & 
Maeder j2000| . We note that all program stars with Mb'i 
magnitudes la 


magnitudes rail well within the ranges in temperature 
and luminosity which are appropriate for RSGs. We cal¬ 
culate stellar masses by interpolating fitted temperatures 
and gravities and calculated luminosities to the Geneva 
evolutionary models. Errors in mass are calculated with 
a monte carlo simulation in which, for 1000 trials, we add 


noise to Teff, log 5 , and luminosity on the scale of our fit 
uncertainties. The mass values and errors presented in 
Table represent the median and Icr standard deviation 
of those monte carlo experiments. This experiment indi¬ 
cates an initial mass range between 12 and 20 Mq. 


5. CONCLUSIONS AND FUTURE WORK 

The quantitative spectroscopy of individual stars in 
distant galaxies as a way to constrain the evolution of 
galaxies has long been a dream of stellar astronomers. 
We believe that with the results presented here we have 
made a major step forward to turn this dream into real¬ 
ity. The stunning agreement between two completely in¬ 
dependent stellar spectroscopic methods analyzing mas¬ 
sive stars in different evolutionary phases with different 
sets of model atmospheres, spectral ranges, ionization 
stages and spectral lines clearly indicates that extra- 
galactic stellar spectroscopy has matured to become an 
accurate and powerful tool. With new multi-object NIR 
spectrographs such as KMOS at the VLT and MOSFIRE 
at Keck and the “classical” MOS instruments at visual 
light for BSGs such as FORS/VLT and LRIS/Keck we 
have now two reliable tools to investigate the chemical 
evolution of galaxies out to ^7 Mpc. With future tele¬ 
scopes with 30 to 40m apertures as the TMT or E-ELT 
and adaptive optics (AO) supported NIR multi-object 
spectrographs t he RSG J-band me thod can be pushed 
much further. Evans et al.| (2011) have demonstrated 
that for AO supported KMOS-like spectrographs magni¬ 
tudes down to J^24 mag could be reached with sufficient 
signal-to-noise in two nights of observing. This will al¬ 
low studies like the one presented here to be carried out 
at 50 Mpc distance and will make a significant volume 
of the local universe accessible for accurate quantitative 
investigations of the chemical evolution of galaxies. 
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TABLE 3 

Comparison of NGC 300 Metallicity and Gradient Literature 


Study Gentral Abundance Metallicity Gradient Notes 

dex, metals 12+log(0/H) R/R25 dex kpc“^ 
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-0.07 ±0.09 

-0.44 ± 0.06 

-0.081 ± 0.011 

Blue Supergiants, Metals 

tsresoim et ai.||Ajuy|l 

8.57 ± 0.02 

-0.41 ± 0.03 

-0.077 ± 0.006 

Hii regions, auroral oxygen 


This work —0.03 ± 0.05 —0.44 ± 0.08 —0.083 ± 0.014 Red Supergiants, Metals 


Note. — Three independent measurements of the evolution of chemical abundance across the star forming disk of NGC 300. 
The agreement between BSGs, RSGs, and auroral Hll region methods is excellent. 
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TABLE 4 


Target 

Alt. Designation 

Right Ascension 
HH:MM:SS 

Declination 

DD:MM:SS 

R/R.25 

HST F814W 
mag 

IRAC Band 1 
mag 

IRAC Band 2 
mag 

010 


00:54:52.07 

-37:40:32.1 

0.07 

18.220 

± 

0.001 

14.61 

± 

0.03 

14.61 

± 

0.05 

022 


00:54:53.83 

-37:40:08.6 

0.12 

18.956 

± 

0.002 

15.88 

± 

0.03 

16.03 

± 

0.05 

024 


00:54:53.12 

-37:39:28.2 

0.20 

18.977 

± 

0.002 

15.78 

± 

0.02 

15.47 

± 

0.07 

007 

100544792-3739105=" 

00:54:47.90 

-37:39:10.0 

0.25 

17.858 

± 

0.001 

15.35 

± 

0.02 

15.45 

± 

0.05 

013 


00:54:49.58 

-37:38:50.5 

0.28 

18.534 

± 

0.001 







127 

100544267-3739172=" 

00:54:42.64 

-37:39:16.9 

0.29 




14.71 

± 

0.02 

14.76 

± 

0.04 

Oil 


00:54:49.91 

-37:38:27.3 

0.33 

18.471 

± 

0.001 

15.90 

± 

0.02 

15.84 

± 

0.04 

009 

MPG2004 72*" 

00:54:45.40 

-37:38:35.2 

0.33 

18.190 

± 

0.001 

15.66 

± 

0.04 

15.70 

± 

0.05 

014 


00:54:48.56 

-37:38:02.0 

0.38 

18.617 

± 

0.001 

15.46 

± 

0.02 

15.29 

± 

0.06 

006 


00:54:47.96 

-37:37:39.2 

0.43 

17.782 

± 

0.001 

15.67 

± 

0.02 

15.72 

± 

0.04 

101 


00:55:02.12 

-37:37:50.3 

0.49 




15.28 

± 

0.02 

15.30 

± 

0.05 

102 


00:54:53.72 

-37:36:47.5 

0.55 




15.22 

± 

0.02 

15.30 

± 

0.04 

133 

100543092-3737570=" 

00:54:30.96 

-37:37:57.0 

0.56 




15.57 

± 

0.02 

15.58 

± 

0.05 

130 

100542809-3738431=" 

00:54:28.10 

-37:38:42.8 

0.57 










126 


00:54:25.54 

-37:39:54.6 

0.59 




15.22 

± 

0.02 

15.26 

± 

0.02 

135 


00:54:29.38 

-37:37:55.9 

0.59 




15.26 

± 

0.02 

15.35 

± 

0.06 

134 

MPG2004 31*" 

00:54:28.25 

-37:37:54.3 

0.61 










039 


00:54:25.76 

-37:37:58.7 

0.65 

18.367 

± 

0.003 

15.99 

± 

0.02 

15.83 

± 

0.07 

128 


00:54:22.00 

-37:39:08.0 

0.67 




14.91 

± 

0.04 

14.87 

± 

0.05 

129 

100541708-3738551=" 

00:54:17.06 

-37:38:54.6 

0.78 










139 


00:54:26.77 

-37:35:04.3 

0.87 

18.283 

± 

0.001 







034 

MPG2004 15*" 

00:54:18.59 

-37:35:38.4 

0.92 

18.624 

± 

0.001 







055 


00:54:22.70 

-37:34:40.0 

0.95 

19.012 

± 

0.001 







132 


00:54:08.74 

-37:38:02.9 

0.96 










140 

100541835-3735024=" 

00:54:18.34 

-37:35:02.0 

0.97 










031 


00:54:23.27 

-37:34:20.3 

0.97 

18.358 

± 

0.001 







037 


00:54:16.18 

-37:34:51.6 

1.01 

18.907 

± 

0.002 








Note. — For objects with no HST F814W entry HST photometry was either not available or uncertain. 
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TABLE 5 


Target Mass Luminosity 
M© log L/L© 


010 

17.4 

± 

1.0 

4.92 

± 

0.08 

022 

12.3 

± 

0.8 

4.64 

± 

0.08 

024 

14.3 

± 

0.8 

4.64 

± 

0.08 

007 

18.8 

± 

1.1 

5.05 

± 

0.08 

013 

16.1 

± 

0.9 

4.80 

± 

0.08 

127 







Oil 

16.3 

± 

0.9 

4.82 

± 

0.08 

009 

17.5 

± 

1.0 

4.93 

± 

0.08 

014 

15.5 

± 

0.9 

4.77 

± 

0.08 

006 

19.2 

± 

1.2 

5.08 

± 

0.08 

101 







102 







133 







130 







126 







135 







134 







039 

16.6 

± 

0.9 

4.86 

± 

0.08 

128 







129 







139 

16.9 

± 

1.0 

4.89 

± 

0.08 

034 

15.8 

± 

0.8 

4.77 

± 

0.08 

055 

14.4 

± 

0.8 

4.62 

± 

0.08 

132 







140 







031 

16.7 

± 

1.0 

4.86 

± 

0.08 

037 

14.8 

± 

0.8 

4.66 

± 

0.08 


Note. 



























